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In the standard cosmological model, the dimming of distant Type la supernovae is explained by 
invoking the existence of repulsive 'dark energy' which is causing the Hubble expansion to accelerate. 
However this may be an artifact of interpreting the data in an (oversimplified) homogeneous model 
universe. In the simplest inhomogeneous model which fits the SNe la Hubble diagram without dark 
energy, we are located close to the centre of a void modelled by a Lemaftre-Tolman-Bondi metric. 
It has been claimed that such models cannot fit the CMB and other cosmological data. This is, 
however, based on the assumption of a scale-free spectrum for the primordial density perturbation. 
An alternative physically motivated form for the spectrum enables a good fit to both SNe la (Con- 
stitution/Union2) and CMB (WMAP 7-yr) data, and to the locally measured Hubble parameter. 
Constraints from baryon acoustic oscillations and primordial nucleosynthesis are also satisfied. 

PACS numbers: 98.80.Es,98.65.Dx,98.62.Sb 



I. INTRODUCTION 

The simplest cosmological model consistent with the 
spatial flatness expectation of inflation is the Einstein- 
de Sitter (EdS) universe with Q m — 1 and Q.k = 0. 
This formed the basis of the "standard cold dark mat- 
ter" (SCDM) cosmology which provided a good descrip- 
tion of the large-scale structure of the universe [3, H[ . It 
was noted however that it is inconsistent with the angular 
power spectrum of the clustering of galaxies in the APM 
survey [3j]. Subsequently large-angle anisotropies in the 
cosmic microwave background (CMB) were detected by 
COBE [4jj thus providing an absolute normalization of 
the amplitude of primordial density perturbations. The 
SCDM model assumes that the power spectrum of the 
primordial density perturbations has the scale-invariant 
form: P{k) = \Sk\ 2 oc fc™, with n — 1 and the predicted 
amplitude of matter fluctuations on small (cluster and 
galaxy) scales is then too high relative to observations 
[5[ . However if the epoch of matter-radiation equality is 
delayed by lowering the CDM density to O m ~ 0.3, then 
the peak in the power spectrum of density fluctuations 
is shifted to larger scales thus decreasing the power on 
small scales and enabling a match to the data. To main- 
tain spatial flatness a non-zero value of the cosmological 
constant was invoked, with A ~ 27Jq corresponding to 
f?A = A/3i?o ~ 0.7. Subsequently it was also observed 
that Type la supernovae (SNe la) at redshift z ~ 0.5 
appear ~ 25% fainter than expected in an EdS universe 
@> 01- Together with measurements of galaxy cluster- 
ing in the 2dF survey || and of cosmic microwave back- 
ground (CMB) anisotropies by WMAP [Sf, this changed 
the standard cosmological model to an accelerating uni- 
verse with a dominant cosmological constant term, which 
has been widely interpreted as a manifestation of the 
physical vacuum or "dark energy" . This "concordance" 
ACDM cosmology (with n A ~ 0.7, Q m ~ 0.3, h ~ 0.7) 
has proved to be consistent with other cosmological data, 
in particular, baryonic acoustic oscillations detected in 
the SDSS 10] and measurements of mass fluctuations 



from clusters and weak lensing [llfl . Further observations 
of both SNe la [jMM and the WMAP 3-year results [H 
have continued to firm up the model. 

Embarrassingly, however, this model lacks a physical 
basis. There are two serious problems with the notion 
that the universe is dominated by some form of vacuum 
energy. The first is the notorious fine-tuning problem 
of vacuum fluctuations in quantum field theory — the 
energy scale of the inferred cosmological energy density is 
~ 10 -12 GeV, which is many orders of magnitude below 
the energy scale of ~ 10 2 GeV of the Standard Model 
of particle physics, not to mention the Planck scale of 
M P = (SttGn)- 1 / 2 ~ 2.4 x 10 18 GeV 0. The second 
is the equally acute coincidence problem: since n\/n m 
evolves as the cube of the cosmic scale factor a(i), there 
is no reason to expect it to be of 0(1) today, yet this is 
supposedly the case. 

Interestingly the WMAP results alone do not require 
dark energy if the assumption of a scale-invariant primor- 
dial power spectrum is relaxed. This is well justified given 
our present ignorance of the physics underlying inflation 
which is believed to have created these fluctuations. It 
has been demonstrated [13, [H| that the temperature an- 
gular power spectrum of an EdS universe with h ~ 0.44 
matches the WMAP data well if the primordial power 
is enhanced by ~ 30% in the region of the second and 
third acoustic peaks (corresponding to spatial scales of 
k ~ 0.01 - O.lh Mpc -1 ). This alternative model with 
no dark energy has a slightly better \ 2 f° r the fit to 
WMAP-3 data than the concordance "power-law ACDM 
model" and, inspite of having more parameters, has an 
equal value of the Akaike information criterion used in 
model selection. Other EdS models with a broken power- 
law spectrum [lj| have also been shown to fit the WMAP 
data. Moreover, an EdS universe can fit measurements of 
the galaxy power spectrum if it includes a ~ 10% compo- 
nent of hot dark matter in the form of massive neutrinos 
of mass ~ 0.5 eV p~7i - FT9j . Clearly the main evidence for 
dark energy comes from the SNe la Hubble diagram. 

It should be kept in mind that the acceleration of the 
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expansion rate is not directly measured but inferred from 
measurements of the apparent magnitudes and redshifts 
of SNe la. Indeed all the evidence for dark energy is 
geometrical, i.e., based on interpreting the data in an 
assumed homogeneous model universe. So far there is 
no convincing observation of dynamical manifestations, 
e.g., the "late integrated Sachs- Wolfe effect". In fact 
what is actually inferred from observations is not an en- 
ergy density, just a value of O(Hq) for the otherwise un- 
constrained A term in the Friedmann equation. It has 
been suggested that this may simply be an artifact of 
interpreting imprecisely measured cosmological data in 
the oversimplified framework of an universe assumed to 
be described by the exactly isotropic and homogeneous 
Friedmann-Robertson- Walker (FRW) metric, in which 
H ~ 10~ 42 GeV - (KFcm)- 1 is the only scale [H. 
Note that the non-zero value of 57a = A/3-ffp is inferred 
from the "cosmic sum rule" VL m + f2/< + = 1, which 
is just a restatement of the Friedmann equation. If how- 
ever this equation does not describe the real universe ex- 
actly, and in order to do so other non-zero terms ought 
to have been added to the sum rule, then we may mis- 
takenly infer a value for J7a of 0(1) if these other terms 
are in fact important. For example, in an inhomogeneous 
universe averaged quantities satisfy modified Friedmann 
equations which contain extra terms since the operations 
of spatial averaging and time evolution do not commute 
2l| . These "backreaction" terms depend upon the vari- 
ance of the local expansion rate and hence increase as 
inhomogenities develop. However although backreaction 
behaves just like a cosmological constant, whether its ex- 
pected magnitude can indeed account for the apparent 
cosmological acceleration is debated and remains an open 
question at present [H-[3l|. 

Another possibility is that inhomogeneities affect light 
propagation on large scales and cause the luminosity 
distance-redshift relation to resemble that expected for 
an accelerating universe. This has been investigated for 
a "Swiss-cheese" universe in which voids modelled by 
patches of Lemaftre-Tolman-Bondi (LTB) space-time are 
distributed throughout a homogenous background. How- 
ever, the results depend on the specific model: some au- 
thors find the change in light pro pag ation to be negligi- 
ble because of cancellation effects [3214351 . whereas others 
claim it can partly mimic dark energy |36H38l | . It may be 
that observers preferentially choose sky regions with un- 
derdense foregrounds when studing distant objects such 
as SNe la, so the expansion rate along the line of sight 
is then greater than average; such a selection effect may 
also allow an inhomogeneous universe to fit the observa- 
tions without dark energy [39| . 

In this paper we are mainly interested in a "local void" 
(sometimes referred to as "Hubble bubble") as an ex- 
planation for dark energy; to prevent an excessive CMB 
dipole moment due to our peculiar velocity we must be 
located near the centre of the void. An underdense void 
expands faster than its surroundings, thus younger su- 
pernovae inside the void would be observed to be reced- 



ing more rapidly than older supernovae outside the void. 
Under the assumption of homogeneity this would lead 
to the mistaken conclusion that the expansion rate of 
the universe is accelerating, although both the void and 
the global universe are actually decelerating. The local 
void scenario has been investigated by several authors us- 
ing a variety of methods [40l - [65j . By modelling the void 
as a open FRW region joined by a singular mass shell 
to a FRW background, it was found that a void with 
radius 200 Mpc can fit the supernova Hubble diagram 
without dark energy [46| . It was also shown that a LTB 
region which reduces to a EdS cosmology with h = 0.51 
at a radius of 1.4 Gpc can match both the supernova 
data and the location of the first acoustic peak in the 
CMB [5l| . Ref.[59| attempted to find the smallest possi- 
ble void consistent with the current supernova results — 
their LTB-based 'minimal void' model has a radius of 350 
Mpc. Unfortunately, since this model is equivalent to an 
EdS universe with h — 0.44 outside the void where the 
SDSS luminous red galaxies lie, as it stands it is unable 
to fit the measurements of the baryonic acoustic oscilla- 
tion (BAO) peak at z ~ 0.35 [6^|. LTB models of much 
larger voids were considered in Ref. [62j (with radii of 2.3 
Gpc and 2.5 Gpc and Hubble contrasts of 0.18 and 0.30 
respectively) and it was demonstrated they can fit the 
BAO data, as well as the SNe la data and the location 
of the first CMB peak. Ref. 64] found the best fit to the 
SNe la data for a void of radius 1.3±0.2 Gpc and Ref.[65| 
confirmed that such a void provides an excellent fit to the 
"Union" dataset of SNe la. 

In this paper we demonstrate that, contrary to the 
results obtained in Refs. [67l - l70l ]. a Gpc-sized void can si- 
multaneously fit the SNe la data as well as the full CMB 
power spectrum, while also satisfying constraints from 
local Hubble measurements, primordial nucleosynthesis 
and the BAO data, if the primordial power spectrum is 
not assumed to be nearly scale-invariant. The layout 
of the paper is as follows. In Sec. [II] we summarize the 
general relativistic framework for LTB models and de- 
scribe the characterization of the void. In Sec. Mil we 
discuss the form of the primordial power, and present a 
physical model with a primordial power spectrum that 
is not scale-free. In order to compare observables in the 
void model to existing cosmological data, some formal- 
ism needs to be developed. This is done in Sec. IIV1 and 
the statistical approach is discussed in Sec. [V] Finally 
Sec. |VT] presents the main results of the paper. 



II. LTB VOID MODELS 

A. The metric and solution 

We model the void as an isotropic, radially inhomoge- 
neous universe described by the LTB metric: 

ds 2 = ~c 2 dt 2 + A ' 2{ ^\ dr 2 + A 2 (r, t)dQ 2 , (1) 
1 + K{r) v ' w 
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where a prime denotes the partial derivative with respect 
to coordinate distance r, and the curvature K(r) is a 
free function, bounded by K < 1. This reduces to the 
usual FRW metric in the limit where A(r, t) — > a(t)r and 
K(r) — > nr 2 . 



We define two Hubble rates: 



rr _A(r,t) 



H 



A'(r,t) 
~ A'(r,ty 



(2) 



where an overdot denotes the partial derivative with re- 
spect to t. The analogue of the Friedmann equation is 



H 2 = 



F(r) c 2 K(r) 



A 3 (r,t) A 2 {r,ty 



(3) 



where F(r) > is another free function which determines 
the local energy density through 



8irGp(r,t) 



F'(r) 



A 2 {r,t)A'{r,t)' 



(4) 



We define dimensionless density parameters Q.m{t) and 
f2if(r) such that 



and 



F(r) = H 2 (r)n M (r)Al(r), 



c 2 K(r)=H 2 (r)n K (r)A 2 (r), 



(5) 



(6) 



where Ho(r) and Ao(r) are the values of H±(r,t) and 
A(r, t) respectively at the present time t = to- The Fried- 
mann equation then becomes [7lj |: 
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(7) 



so £Im ( r ) + ^K{ r ) — 1- This equation can be integrated 
from the time of the Big Bang, fa = ts(r), to yield the 
age of the universe at any given (r, t): 



t - t B (r) 



A/A 



dx 



Ho(r) Jo ^/n M (r)x-^ +tl K {r) 



(8) 



We thus have two functional degrees of freedom, in VLm (t) 
and t&{r), which can be chosen as desired. (The third 
function, ^4o( r )j corresponds to a gauge mode and we 
choose to set Ao(r) — r.) A spatially varying t-Q corre- 
sponds to a decaying mode [72j |. so for simplicity we set 
<b = everywhere, so that at the current time to'- 



Ho(r) 



-y/—Q K + £l M sin 



to (-Sl K ) 
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n K - fi M sinh _1 JS*- 
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, fix < 0; 
= 0;(9) 

fiif > 0. 



The void model can then be specified by the choice of one 
free function, which we take to be f2jf(r), and a constant 
H = Ho(0) which determines the local Hubble rate at the 
centre (and is equivalent to choosing to). Note that both 
at the centre of the void and far outside the void, the 
definition of fl m (r) reduces to the standard FRW density 
parameter f2 m . 

The solution to Eq.Q for general r and t can be given 
in parametric form for the different values of J7j<-(r) (or 
K{r)) as follows 0: 



• for Cl K (r) > 0: 

n M (r)A (r 



A 



Hot = 



2Q K (r) 



(coshr/-l), (10a) 
(sinh^-ry). (10b) 



for tt K {r)= 0: 

A=\ (18Q M (r)) 1/3 (H tf 3 A (r). (11) 



for n K (r) < 0: 



A= ^M(r)A (r) (1 _ cosu) 



Hot 



2\n K (r)\ 
2\n K (r)\V 2 



(u — sinu) 



(12a) 
(12b) 



Light travels to an observer at the centre of the void 
along null radial incoming geodesies described by [42j |: 

dt 



1 



dz {l + z)H\\(zy 
dr cy/l + K (r) 



dz (l + z)A'(z)H l \(z)' 



(13a) 



(13b) 



where Hu(z) = Hu (r(z),t(z)) etc. The angular diameter 
distance at redshift z is then given by 



d A (z) = A(r(z),t(z)), 
and the luminosity distance by 

d h (z) = (l + z) 2 A(r(z),t(z)). 



(14) 



(15) 



All observable quantities along the light cone can be cal- 
culated from these equations. 



B. Void profile 

We can now choose the void profile by specifying 
n^(r) and the local value of the Hubble rate by speci- 
fying H. Although the void profile may have any shape, 
we restrict ourselves to the simple Gaussian form: 



fijtfW = ^out - (fiout - ^in) exp 



(16) 
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where fii n and f2 ut correspond to the matter density pa- 
rameter at the centre of the void and at infinity, respec- 
tively, and ro characterizes the width of the void. We 
wish to look only at voids that are asymptotically EdS, 
so we set f2 ut = 1 ■ Thus distances in the void model are 
completely specified by the three parameters fii n , ro and 
H. 

We wish to stress that in restricting ourselves to voids 
which have a Gaussian profile, we may be missing the 
model that fits the data best. In principle we could sam- 
ple a wider class of profiles and choose the form that 
gives the best fit. However in this paper we are concerned 
mainly with providing a counter-example of a void which 
can simultaneously fit both the SNe la magnitudes and 
the CMB spectrum, so we do not perform this search. 



III. PRIMORDIAL POWER SPECTRA 

The observed power spectrum of CMB anisotropics 
in any cosmological model is a convolution of three un- 
knowns. The first two are the local physics at the time 
of recombination, which is dependent on the composition 
of the universe at that time, and the angular diameter 
distance to the last scattering surface (LSS), which de- 
pends on the geometry of the universe. Both of these are 
completely specified by the choice of the void model as 
described in Section fll B[ along with the further specifica- 
tion of the baryon fraction, fib and the baryon-to-photon 
ratio 77 = n^/rij. 

The third unknown is the shape of the primordial 
power spectrum of density perturbations. In toy mod- 
els of inflation this is close to scale-invariant and feature- 
less. Under the assumption that the primordial spectrum 
is described by a simple power law, the standard ACDM 
concordance cosmological model fits the observed angular 
power spectrum reasonably well 73]. However, there is 
no independent evidence for this form of the primordial 
power and as the observed anisotropies arise as a con- 
volution of the assumed primordial spectrum with the 
transfer function of the assumed cosmological model, it 
is clear that we cannot determine one without making 
assumptions about the other. (In fact, there are indica- 
tions that the primordial spectrum is not scale-free, even 
when a ACDM cosmology is assumed [t3 - [83| .) 

In Refs. (67l - [7Cj| it is argued that, assuming near-scale 
invariance of the primordial power, void models cannot 
simultaneously provide an explanation for supernovae 
magnitudes and fit the observed CMB spectrum, unless 
the void is extremely deep [6jJ or embedded in a universe 
with large non-zero overall curvature [70l ]. However, in 
considering void models as an alternative to dark energy, 
we are in any case departing from the concordance cos- 
mology, so there is no need to retain the assumption that 
the primordial power spectrum is scale-free. If this as- 
sumption is relaxed then even an EdS cosmology without 
dark energy can fit the CMB data [Tt], Ell . In this paper 
we consider a particular alternative form of the primor- 



dial V(k) that is not scale- free but is in fact physically 
well-motivated, as described below. 

A. Bump model 

Whereas the simplest toy models of inflation contain 
only a single scalar field which rolls slowly down its po- 
tential, physical models gcncrically contain other fields, 
whose evolution is typically not slow-roll. These fields 
may couple to the inflaton and affect its evolution, thus 
breaking the scale-free nature of the primordial power, 
with important consequences. 

An ex amp le of such a physical model is "multiple in- 
flation" [84[ in the framework of N = 1 supergravity, 
the locally realised version of supersymmetry (SUSY). 
This model includes "flat-direction" fields ip which have 
gauge and Yukawa couplings to ordinary matter but are 
only gravitationally coupled to the inflaton. Such flat- 
directions have been classified and tabulated in the Min- 
imal Supersymmetric Standard Model (MSSM) [Hj]. As 
the universe cools during inflation, the ip fields undergo 
symmetry-breaking phase transitions, causing a sudden 
change in the effective mass of the inflaton (which is as- 
sumed to be a field in a "hidden sector"). For a single 
flat-direction field, this produces a 'step' in the spectrum 
of the curvature perturbation (86j . If more than one flat- 
direction field is present, they can couple to the infla- 
ton with opposite signs, and instead produce a 'bump' 
feature in the power spectrum [171 ] . (A toy model that 
produces a similar "step" feature in the power spectrum 
has also been proposed [13] and has been studied with 
respect to fitting the WMAP 1-year [H, HU and 3-year 
90] data, albeit in a different context to that considered 
here. Other signatures of multiple inflation, in particular 
the generation of associated non-Gaussianities, have also 
been studied (9l| .) 

The potential for the inflaton <fi and the flat-direction 
fields ipi is then given by 

V{</>, V>i> 2 ) =V - \m 2 ttU 2 + X -\ X HU 2 ^\ 

+ 7i^i 1 + 

-\t4Hfy 2 2 +j2^, (17) 

where mH\ and fiiHj are the masses of the <f> and ipi fields 
respectively, XiHf is the coupling of the ipi field to the 
inflaton, 7, is the coefficient of the non-renormalizable 
operator of order n, which lifts the potential of the ipi 
field and H\ is the Hubble scale during inflation (units 
of Planck mass Mp = 1 are used throughout). The two 
flat-direction fields remain trapped at the origin by ther- 
mal effects until the phase transitions take place at times 
tj. Natural values for the parameters of this model are 
I Ai| ~ 1, /if ~ 3, and 7, ~ 1. According to the list of 
flat-direction fields in Ref. (85|, the ones with rij = 12, 16 
will be the most relevant. It is usually assumed that 
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some symmetry protects the inflaton mass m from SUSY- 
breaking corrections (the "77 problem"), thus allowing 
sufficient e-folds of inflation to occur. Note that the times 
of the phase transitions t\ and t% — or, equivalently, the 
scales k\ and &2 at which the effects of the flat-direction 
fields begin to be felt — are however arbitrary. 

In Ref. [13] a full likelihood analysis was performed to 
demonstrate that a primordial power spectrum with such 
a 'bump' would allow even an EdS universe with A = 
to fit the WMAP 3-year data, although this requires a 
rather low value of ho ~ 0.44. As void models generically 
imply a lower global value of ho than the local value at 
the centre where we are located, such a primordial power 
spectrum could help a void model to fit the CMB. 

The primordial power spectrum produced by this 
'bump' model can be calculated exactly by numerically 
solving the equation of motion of the cf> field according to 
the potential in Eq. (fTT|) . However, we do not choose to 
do so for two reasons. Firstly, this evaluation is computa- 
tionally expensive and slows an MCMC analysis consid- 
erably. Secondly, our results regarding void models are 
dependent only on the existence of such a feature, rather 
than the details of the particular theory that explains it. 
Multiple inflation provides an example of such a theory, 
but similar features with steps and oscillations may also 
be produced, e.g., in DBI inflation (92j. Strictly speak- 
ing it is not even necessary to assume that these features 
arise from an inflationary model at all. 

With this in mind, we introduce a parameterization of 
the primordial scalar power V(k) to capture the essential 
features of the 'bump' model instead of attempting to 
reproduce it exactly: 



V(k) = Vo ( 1 - atanh(6a;) + ce" 



-{bxy 



(18) 



where x = \og 10 (k/ko). In Fig. Q]we plot V(k) obtained 
from a full calculation in the "bump" model, with param- 
eter values ni = 16, n% = 12, [i\ = ju| = 3, 71 = 72 = 1, 
Ai = —0.2, A 2 — 1 and m 2 = 0.05 (all set at their "nat- 
ural" values, save for Ai). Hi, k\ and ki are free pa- 
rameters and are chosen so as to broadly reproduce the 
best-fit primordial spectrum from Ref. [17] . In the same 
figure we also plot the form of our simple parameteriza- 
tion. We do not consider any tensor modes, as in "small 
field" models such as multiple inflation these are always 
negligibly small. 

This form of V{k) is thus determined by 5 free param- 
eters: Vo, a, b, c and ko- This may appear to be a step 
backwards from the standard power-law form, which de- 
pends on only 3 (As, n s and n r un — the amplitude, slope 
and "running" of the spectrum) . However it is important 
to bear in mind that we are merely using an empirical pa- 
rameterization for a V(k) that is ultimately the result of 
some underlying theory, in which all of these parameters 
are in fact determined by fundamental physics. 

On the other hand, it can be argued that the standard 
parameterization of a power-law hides other parameters 
which are artificially set to zero through the overriding 
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FIG. 1. The solid (blue) curve shows the primordial power 
spectrum obtained from a full calculation of the multiple 
inflation model, with parameters as described in the text. 
The dashed (red) curve shows the simple parameterization 
(Tl8)) . with parameter values T> = 2.48 x 10~ 9 , a = 0.1, 
b = 2.0, c = 3.0 and ko = 0.015. For comparison, the dot- 
ted (black) curve shows a standard power-law spectrum with 
slope n B — 0.954, amplitude As = 2.3 x 10 -9 , which are the 
best-fit values for the ACDM model with no "running" to the 
CMB+SN cf A datasets (Sec. ED- 



assumption that only a single scalar field is involved. In 
the absence of an accepted physical mechanism for infla- 
tion, we feel that all plausible alternatives ought to be 
considered, rather than judging purely on the basis of 
naive parameter-counting. It is with this rationale that 
we use this form of the primordial power in constructing 
our counter-example. 



IV. FITTING THE MODEL TO OBSERVATIONS 

In this section we will discuss how to fit the void model 
to the best available cosmological data: the WMAP 7- 
year release [73| , SNe la data from the Constitution [§3| 
and the Union2 compilation (94[, the local Hubble rate 
measurement by the Hubble Key Project [95| , constraints 
from big bang nucleosynthesis [96| . and BAO measure- 
ments from SDSS DR7 [<tfj]. 

As cosmology in an LTB metric differs from the stan- 
dard FRW approach, we need to develop some formal- 
ism to allow us to confront these datasets in a consistent 
manner. This is outlined in the following. 



A. SNe la magnitudes 

The distance modulus of a supernova is defined as the 
residual between its apparent and absolute magnitudes 
m and M, and is related to its luminosity distance as 



fi = m — M = 5 logj 



Mpc 



25, 



(19) 
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so that, knowing fi for each object in a dataset, we can 
use Eq. (fT5)l to compare a void model to the data. 

We use two datasets: the Constitution compilation 
of 397 SNe in the redshift range z = 0.015 - 1.55 [H 
and the Union2 compilation of 557 SNe in the range 
z = 0.015 — 1.4 [94j . The Constitution sample uses 
the SALT lightcurve fitter while Union2 uses the newer 
SALT2 fitter. Neither fitter, however, directly provides 
the value of \i for each object — for instance, the SALT 
fitter provides as output values of m^ ax (the rest-frame 
peak magnitude in the B-band), a time stretch factor, 
s, and a colour parameter, c, for each supernova, from 
which the distance modulus is calculated as: 



m^T - M + a ■ { Si 



l)-0 



(20) 



where M, a and /3 are empirical coefficients whose values 
are determined by marginalizing over the fit to a partic- 
ular fiducial cosmology. The SALT2 fitter follows the 
same principle, though with a 'stretch factor', x%, that is 
analogous to but not the same as s. 

Both the Constitution and Union2 compilations pro- 
vide tabulated values of distance modulus fj, and its es- 
timated error oy, for a choice of coefficients M , a and (3 
calculated for the best-fit flat ACDM model. In addition, 
includes an important contribution from a systematic 
uncertainty whose value is chosen such that the reduced 
X 2 value for the flat ACDM model is of order unity. Hence 
when comparing any alternative model to ACDM on the 
basis of their respective fits to the SNe la data, we ought 
to in principle perform the entire fitting analysis from 
scratch, choosing the values of M, a and /3 that produce 
the best fit for the particular model under consideration, 
and adjusting the systematic error inserted by hand in an 
equivalent manner. However, to do so we would require 
not only the actual lightcurve fitter outputs for each su- 
pernova (which the publicly downloadable data for the 
Constitution set does contain, but Union2 does not), but 
also the full covariance matrix of the lightcurve fits, which 
neither compilation provides^ Therefore, we are forced 
to adopt the common (if incorrect) procedure of fitting 
the void model by marginalizing over the unknown ab- 
solute magnitude M alone, while bearing in mind the 
caveat that this may bias our results and our best-fit 
void parameters in an unknown manner. 



B. CMB power spectrum 

It has been shown in several previous studies that void 
models are only very loosely constrained by the posi- 
tion of the first peak of the CMB [51, 62]. However, the 
WMAP satellite has in fact measured the angular power 
spectrum Qs over a wide range of multipoles, I, and it is 



1 This is provided in the SNLS 3-year data release [iisll ; however it 
contains only 231 objects, so we do not consider it here. 



preferable to use as much of this data as possible. This 
entails the calculation of the full predicted Q spectrum 
for the void model. This has been done using various 
methods f67j - l7"ol ] . but always with the assumption of a 
nearly scale-invariant primordial spectrum. 

In order to be able to calculate the predicted angu- 
lar power spectrum for an LTB model using one of the 
publically available Boltzmann codes (we use a version 
of C AMB 99] , modified to accept the primordial power 
spectra we consider), we use a version of the "effective 
EdS approach" (68, 69]. In brief, this consists of con- 
structing an effective EdS model which has the same 
physics at recombination as the given void model, and 
the same angular diameter distance to the LSS. Then, 
given the same primordial power spectrum, the C;s of 
the effective model and the void model will be the same, 
at intermediate and small angular scales. Thus the effec- 
tive model can be used as a calculational tool to obtain 
the power spectrum for the void model under consider- 
ation. Note however that the effective EdS model will 
in general have a central temperature, TJ^ ds , and a cen- 
tral Hubble rate, -ff EdS , that are different from the actual 
T and H that we observe today. This is because the 
physics has been matched at early times; at late times 
the models must then necessarily differ. 

To generate the effective model, we adopt the follow- 
ing procedure. First we specify the void profile and local 
Hubble rate by choosing the parameters f2; n , tq and H. 
Then, using Eqs. I13ap and (I13b[) we numerically inte- 
grate out from the centre of the void along the past light 
cone to obtain the coordinates {jm , t m ) at an interme- 
diate redshift z m . As in Refs. [68|, [69| we choose this 
redshift to be z m = 100, where the spatial curvature of 
the void is negligible but the radiation density is still 
small and can justifiably be ignored in the calculations. 

At these coordinates we now calculate the Hubble rate, 
H m = H±/u (r m ,t m ) and the angular diameter distance 
A m = A (r m ,t m ). Let us denote by r EdS (t) the comoving 
radial coordinate of a radial light ray in the EdS universe 
and by z Eds the redshifts seen by an observer at r Eds = 0. 
To ensure the matching of the distances to the LSS, we 
impose the condition A m = a(?" Eds )r Eds for EdS scale 
factor a. This provides us with the relation 



An 



where 



2c 



1 



(l + zEd S) tfEdS 



0" 



,EdS 



(21) 



y EdS 



is the redshift in the EdS universe at co- 
ordinates (r Eds ,t m ). (Note that our procedure implies 



y EdS 



7^ z m and thus differs slightly from the equivalent 
method outlined in Refs. fo8l l69| where instead z EdS =/= 0.) 
To ensure the same physics at early times in the two mod- 
els we also match the Hubble rates, H m = H Eds (z EdS ) , 
and this combined with Eq. (|2ip provides us with an ex- 
pression for z EdS 



,EdS 



4c 2 



(22) 
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Given the value of z Eds , we can then calculate the effec- 
tive EdS mean temperature and Hubble rate via 



(23a) 
(23b) 



rriEdS 
1 


= T ( 






crEdS 
H 






(1- 





where T = 2.726 ± 0.001 K is the CMB temperature 
observed today [100[ . 

The parameters T Ed ^ # Eds , 0„ = l,fi A = Q K = 
and the baryon fraction £lb can then be fed into CAMB, 
together with the choice of primordial power spectrum, 
to find the spectrum of Qs for the effective EdS model. 
This will be identical to the spectrum actually observed 
at the centre of the void, at large enough values of I. 

At small values of / however, the void model will cre- 
ate an integrated Sachs- Wolfe (ISW) signal which is not 
captured by the effective EdS model. However, no rigor- 
ous calculation has yet been made of this expected ISW 
signal. In addition, the void may also in principle have 
a different reionization history than that of the effective 
EdS model. These two effects mean that at small I the Q 
values calculated using the EdS model will differ, in both 
the TT and TE power spectra, from the actual spectrum 
due to the void. 

In order to account for this, we choose to apply a cutoff 
at I = 32 in both the WMAP-7 TT and TE power spec- 
tra. The value Z = 32 is chosen also to coincide with the 
switch-over point for the TT power spectrum at which 
the WMAP likelihood routine switches between the low- 
l, gibbs-sampling likelihood estimation technique and the 
master code for high-Z (for the TE spectrum the switch- 
over point between pixel-based analysis and the master 
code is close by at I = 24). This also allows a more di- 
rect interpretation of the likelihood L in terms of a x 2 
value. We have checked that increasing the cutoff point 
does not materially affect our results. 

Given that in Ref. [13] it is found that an EdS model 
can fit the WMAP data with a primordial power spec- 
trum similar to the 'bump' model we use here, we 
expect to find a good fit to the CMB with ft EdS = 
J ff^ dS /(100kms" 1 Mpc~ 1 ) ~ 0.44 and T Eds ~ T . This 
should be contrasted with the values of /i Eds ~ 0.51 and 
T Eds ~ 3.4 K required with a power-law primordial spec- 
trum mm. 



C. Local Hubble rates 

It has been claimed in previous studies [o^l - [7"(]| that the 
local Hubble rate of void models that fit the CMB and 
SNe la data simultaneously must be very low (as low as 
ho — 0.45 in Refs. [gl, HU) and that this argues against 
void models, since the measured lo cal v alues (at z < 
0.1) are significantly higher (see Ref. jl0l| for a review). 
This can provide an important discriminant against void 
models. Therefore in performing the MCMC analysis, we 



also fit the local Hubble rate for the void model /ig TB = 
if/100 to the Hubble Key Project (HKP) value h = 
0.72 ± 0.08 [95[. We note that there is some variation 
in the value of ho obtained by different groups, ra ngin g 
from /i = 0.623 ± 0.06 102 to h = 0.742 ± 0.036 [lo3 |. 
and the HKP value lies in between these two. 

The SNe la compilations, WMAP-7 observations and 
the HKP value for ho form the primary datasets that we 
use to constrain the void model. We also discuss the fit 
to some other cosmological data below. 



D. Big bang nucleosynthesis 

Our theoretical understanding of the physics of the 
epoch of big bang nucleosynthesis allows us to use ob- 
servations of the abundances of various elements to con- 
strain the baryon-to-photon ratio, rj — n^/n^, at that 
time. Although r\ is constant with time in FRW space- 
time, this is not the case in LTB spacetime. However, 
assuming that 77 is spatially constant in the LTB mod el 
(which need not necessarily be true, see Refs. |104l . ll05j ). 
we can use the effective EdS model model to calculate 77 
for the void model. Since both models share the same 
early universe physics by construction, we have [69[ 



7710 = 10 10 ?; = 273.9 



To V 

T Ed S J 



(24) 



where Uh = (ho dS ) ■ The inferred primordial abun- 
dance of deuterium, together with that of helium, pro- 
vides the constraint 5.1 < ?/io < 6.5 at 95% C.L. [96| |. 



E. BAO scale 

The BAO data provided by the SDSS collaboration [§3] 
are essentially measurements of a feature in the correla- 
tion function of the observed galaxy distribution which is 
related to the physical sound horizon at the CMB scale, 
evolved down to the redshift at which the measurement 
is made. The full theory for how these perturbations 
should evolve in LTB spacetimes is unknown and a dif- 
ficult problem, as the background curvature enters the 
Bardeen equation for the gravitational potential. As this 
may vary significantly over scales of ~ 150 Mpc, it can 
potentially add a sig nific ant and as yet unknown distor- 
tion to BAO scales [106| . 

Nevertheless, some efforts have been made to compare 
void models to BAO data under the assumption that the 
evolution of perturbations does not depend on scale (as is 
the case in FRW spacetime) [68H7Cl | . Under this assump- 
tion the difference between LTB and FRW spacetimes 
is simply that as H\\ 7^ H± for r > 0, physical length 
scales in an LTB spacetime evolve differently in radial 
and transverse directions at late times, whereas in FRW 
models the evolution is isotropic. It is not clear that this 
is necessarily a valid assumption to make. However, in 
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order to provide a comparison with the previous stud- 
ies we follow the same prescription as in Ref. [6!|, while 
bearing in mind that a better calculation may lead to 
significant changes in the results. 

The first step is to construct another effective EdS 
model, referred to as the "BAO model", which shares 
the same early universe physics as the void model un- 
der consideration, but unlike in the previous case, need 
not match the angular diameter distance to the LSS. We 
choose in this case to match the central temperatures, 
Tq ao = To, which then provides us with the relation: 



rrBAO 
^0 



(1 



v3/2' 



(25) 



as unlike in the previous Section, we will now have 
z m L ° = z m (where, as before, we choose z m = 100). We 
also match the values of fib and Vt m = 1 at this redshift 
to ensure the same early universe physics. 

As the BAO model and the LTB void model 
have, by definition, the same physics at coordinates 
(r (z m ) , t (z m )) = (r TO , t m ), and since at this redshift the 
LTB spacetime is sufficiently close to FRW, we may con- 
clude that the physical sound horizon in the void model 
at time t m is essentially isotropic and homogeneous, and 
equal to the sound horizon s p (z m ) in the BAO model, 
which may be calculated according to [l07j | : 



(1 + 2) vi + io(«b^) 3/4 



(26) 



In order to evaluate the BAO scales observed at redshift z 
by an observer at the centre of the void, we need to evolve 
this physical sound horizon scale down to coordinates 
(r (z) , t (z)) in the LTB background. Thus the radial 
and transverse physical BAO scales at redshift z will be, 
respectively, 

A' (r (*),*(*)) 



i BAO 



(z) = Sp (z m ) 



A' (r(z),t m ) 
A(r(z) ,t(z)) 



(27) 



(28) 



A(r(z),t m ) ' 

and these can in turn be rewritten in terms of the corre- 
sponding redshift and angular intervals: 

(29) 
(30) 



Az{z) = (l + z)l^ (z)if|| (r(z),t(z)) 
A9(z) = 



BAO 



A(r(z),t(z)Y 

It is these values of Az (z) and AO (z) that are directly 
measurable. As mentioned in Ref. [69|], the redshift scale 
Az (z) is expected to be a stronger discriminator of void 
models than A9(z). However, wh ile e s tima tes of the 
radial BAO scale have been made jl08l Il09j |. the sta- 
tistical si gnificance of these claims has been questioned 
[llOl llllj . Given this ambiguity, we choose not to use 
this data; instead we use constraints on the ratio, 



@{z) 



(^rec) 



(31) 



at redshifts z — 0.2 and z = 0.35 from Ref. [9jJ, where 
z roc is the redshift at recombination and Dy (z) is an 
isotropized distance measure. This is given in a FRW 
model by [lfjj 



nFEW 



zd\ (z) 



H(z) 



1/3 



(32) 



where cLa (z) and H (z) represent, respectively, the angu- 
lar diameter distance and the Hubble rate in the FRW 
model. 

As the BAO measurements are quoted in terms of the 
ratio O (2), it is necessary to convert them into a form 
that can be related to Az (z) and AO (z). It can be shown 
Wdi that 



'(«) = (A6 2 Az) 1/3 = zy 3 e(z). 



(33) 



and this is the measure we use to compare the void model 
to the data. 



V. METHOD 

We perform a likelihood analysis using COSMOMC 
[j~L2| to generate Markov-Chain-Monte-Carlo (MCMC) 
chains to estimate confidence limits on the parameters 
in fitting the model to the data. For each void model, 
specified by fl m , ro and H, we first calculate the effective 
parameters i^ ds and T EdS as described in Sec. lIVBI Wc 
fix the optical depth to the LSS to the WMAP-7 value, 
although this choice is immaterial as our cutoff at I — 32 
means that the data do not constrain reionization in any 
case. Therefore the full set of parameters for the void 
model is: local matter density f2j n , void radius ro, lo- 
cal Hubble value H, baryon fraction Q^, and the power 
spectrum parameters Vo, a, 6, c and fc . (The cold dark 
matter density Q c is set equal to 1 — Ob since we are not 
considering large-scale structure formation here so do not 
need to invoke a small component of hot dark matter to 
match the SDSS power spectrum of galaxy cluste ring as 
in Ref. 17].) These are fed as inputs to CAMB [Itf to 
generate the CMB spectrum for the model. 

Once the output Q values have been obtained from 
CAMB, we can fit the model to the WMAP-7 data [73|. 
the SNe la data (the Constitution [93| and Union2 |94l ] 
datasets are fitted separately), the HKP value for the 
local Hubble rate the BAO data [13 and the BBN 
constraint [961 ] . 

In order to compare the goodness of fit, we perform 
exactly the same fitting procedure for a vanilla ACDM 
model whose parameters are: baryon density f2b^ai,cold 
dark matter density fl c h^, local Hubble rate 100ft-oP an d 



2 The COSMOMC code actually uses the ratio of the sound hori- 
zon to the angular diameter distance as the input parameter, but 
this is equivalent to using ho- 
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FIG. 2. The TT and TE power spectra of the models that best fit the CMB and Consti tution SNela data. The solid (red) 
curve is for the void model, calculated using the effective EdS approach described in Sec. IIVB1 and the dashed (black) curve 
is for ACDM. Binned WMAP-7 data is also shown. 
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FIG. 3. Constraints on the effective EdS parameters for the void model with a bump, from CMB and SNe la data. The solid 
(red) contours show the 1 and 2<r likelihood confidence intervals for the Constitution SNe la data set and the dashed (blue) 
contours are for Union2. 



the power spectrum amplitude, A$, and slope, n s . We 
characterize the best-fit likelihood of the void model by 
the value Ax 2 — —2 ln(£ vo i(j/£ACDM)> which means that 
negative values of A\ 2 favour the void. In Fig. [5] the TT 
and TE power spectra are shown for two sample best- 
fit models. The quantitative results of the analysis are 
discussed in the next section. 



VI. RESULTS 

In this section we shall use SNcfA to denote the Con- 
stitution data set and SNtj2 to denote Union2; CMB 
denotes the WMAP-7 data and HKP the Hubble Key 
Project value for the Hubble rate. 

In Figure [3] we plot the 2D likelihoods for the ef- 
fective EdS model parameters T EdS , F EdS and w b ob- 
tained from the MCMC chains for both CMB+SN C fA 
and CMB+SN U2 data. Clearly, T EdS ~ T and /i EdS ~ 
0.45, which is exactly as we expect, given the results of 



10 




3 4 5 13.5 14 14.5 15 

r (Gpc) Age(Gyr) 



FIG. 4. Marginalized ID likelihoods for the void model with a bump, for the fit to CMB and SNe la data. The solid (red) 
curves are for the Constitution SNe la data set and the dashed (blue) curves are for Union2. 



Ref. [171 . Figure H] shows the marginalized likelihoods of 
the void parameters and the resultant age of the universe, 
for the same datasets. 

For both CMB+SN C fA and CMB+SN U2 we obtain a 
local Hubble value ho ~ 0.60 ± 0.02, which, although 
slightly low, is still consistent with the HKP value of 
ho = 0.72 ± 0.08 to within 2a. It is also much higher 
than the value of ho obtained in previous studies [(38|, 
l69j . This is because in order to fit the CMB with a 
power-law primordial spectrum a value of T Eds ~ 3.4 K 
is required and in order to generate this, the void must 
be surrounded by a large over-dense shell. The presence 
of this shell means the asymptotic value Ho(r ^ ro) ~ 
H^ ds ~ H (0) and hence in order to fit the CMB, H o (0) 
must be small. In contrast, our void profile does not have 
an over-dense shell, and as shown in Figure^ T EdS ~ T . 
This allows us to generate the difference between ho — 
0.60 and h$ dS ~ 0.45. 

Our void model also does not suffer from the "old age" 
problem referred to in Ref. [69]: we find the age of the 
void universe to be 14.4 ± 0.3 Gyr for the CMB+SN C fA 
chains and 14.5 ± 0.3 Gyr for CMB+SN U2 , which are 
significantly less than the values of ~ 18.8 Gyr quoted 
in [69]. This is unsurprising as the age to is related to 
the value of ho by Eq. [51 and reasonable values for one 
ensure reasonable values for the other. 

To enable a comparison of the quality of the fits, we 
show in Table U the A\ 2 values for the best-fit model 
with a void and a spectral "bump" , relative to the stan- 
dard ACDM model, for different combinations of the pri- 
mary constraining data sets. It is clear that the A\ 2 
values in all cases are small given the number of degrees 
of freedom involved, which shows that the void model is 



TABLE I. Best-fit likelihood values for ACDM and the rela- 
tive A\ 2 values for the void model with a spectral 'bump for 
different choices of the fundamental constraining data sets. 
Here SNcfA refers to the Constitution sample and SNu2 to 
Union2. CMB refers to WMAP-7 data and HKP to the Hub- 
ble Key Project value for ho- Constitution data favour the 
void slightly and Union2 data favour ACDM slightly, but nei- 
ther difference is significant given the number of degrees of 
freedom. 



Datasets 


# d.o.f. 


— 2ui(£acdm) 


Ax 2 


CMB 


1936 


5785.0 


+0.9 


CMB + SNcfA 


2333 


6250.4 


-3.3 


CMB + SN U2 


2493 


6315.7 


+ 1.4 


CMB + SN cfA + HKP 


2334 


6250.5 


-0.6 


CMB + SN U2 + HKP 


2494 


6315.9 


+3.4 



perfectly compatible with these datasets. A more quan- 
titative comparison cannot be made because as noted in 
Sec. IIV Al the SNe la data are adjusted to fit a fiducial 
ACDM model and the error bars are tuned to give a x 2 
per degree of freedom of order unity for this model. 

It should be noted that the addition of the HKP con- 
straint increases the Ax 2 value by ~ 2 for Union2 and 
~ 2.7 for Constitution, which reflects the fact that CMB 
and SNe la data favour a value of ho that is consistent 
with, but lower than, the HKP value. Even with this 
increase however, the A% 2 values are still small (and 
in fact the void provides a marginally better fit to the 
CMB+SNcfA data than does ACDM). It is possible that 
a thorough exploration of possible void profiles or use of 
unbiased SNe la data will make this fit better, but as this 
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does not detract from the main results presented here, we 
do not investigate it in this paper. 

A more interesting constraint comes from Big Bang nu- 
cleosynthesis. From our MCMC chains for CMB+SN C fA 
and CMB+SN U2 we find in both cases 7710 = 4.6 ± 0.2, 
which is ~ 2a below the best-fit value. This arises pri- 
marily because u>b is on the low side, as seen in Figure [3] 
When the BBN constraint is added to the MCMC chains, 
we find an increase in A^ 2 of ~ 2.5 relative to the val- 
ues for CMB+SNcfA+HKP and an increase of ~ 3.1 for 
CMB+SN U2 +HKP. It should be noted however that the 
inferred primordial abundance of lithium indicates a sig- 
nificantly lower value of r\ than does deuterium [96] , hence 
a better understanding of the chemical evolution of these 
fragile elements is required before the significance of the 
marginal discrepancy above can be assessed. 

Finally, we find that adding the BAO data results 
in an increase in A^ 2 of 4.9 relative to the value for 
CMB+SN U2 +HKP+BBN and 5.5 relative to the value 
for CMB+SNcfA+HKP+BBN, for an additional two de- 
grees of freedom. The main reason for this increase is that 
the BAO data require a shallower void with Sl; n ~ 0.17, 
whereas the CMB+SN data favour a slightly deeper void 
with Oi n ~ 0.13. Most of the increase in \ 2 arising from 
including the BAO data comes from the poorer fit to 
SNe la data that results. Given the uncertainty in in- 
terpreting xIn that we have already mentioned, as well 
as the assumptions that have gone into our analysis in 
Sec. IIVE1 we cannot draw any definite conclusions from 
these results. However it would appear that the BAO 
data are certainly not inconsistent with a void model. 



nates the distortion, this can be written as [6J 



y 
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az 



(34) 



where r is the optical depth, (3 (z) is the dipole temper- 
ature anisotropy in the CMB observed at redshift z, and 
the integral is taken up to the redshift of reionization z rc . 

The FIRAS instrument on COBE provides an upper 
bound y < 1.5 X 10" 5 (at 2a) [Tl3|. While this can 
constrain some void profiles, voids lacking an overdense 
outer shell are found not to be in tension with this upper 
bound 69]. The Gaussian void profile that we consider 
does not have a prominent overdense outer shell. Since 
we do not expect an interesting constraint, we chose not 
to evaluate the y-distortion, which would be uncertain in 
any case since the exact reionization history (and thus 
z rc ) of the void model is unknown. 

Another direct test of a loc al void is via the kinetic 
Sunyaev-Zel'dovich effect 1_14[. This has been applied 



to cluster kSZ observations (e.g., Ref . [63l. |115{ ) and to 
full-sky data Ill6l from the Atacama Cosmology Tele- 
scope (ACT) jll7| . and is found to place constraints 
on large voids. Allowing an inhomogeneous bang-time 
or correct ly accounting for the effects of radiation (see 
[l04l Il05j ) may potentially weaken these constraints 
while preserving the fits described here. Such an analy- 
sis is beyond the scope of the current paper but will be 
investigated in the future. 

In summary, in this paper we have presented a local 
void model which fits SNe la and CMB data, local Hq 
values, nucleosynthesis constraints and BAO without re- 
quiring dark energy and thus provides a counterexample 
to the claim that dark energy is necessary to fit these 
observations. 



VII. DISCUSSION 
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under the assumption that the dipole anisotropy domi- 
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